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1.  Introduction 


Auroral  oval  boundary  has  been  eonsidered  a  eritieal  requirement  for  the  DMSP  mission  and 
spaee  situation  awareness.  The  auroral  boundary  speeifieation  supports  requirements  for  global 
situational  awareness  of  geomagnetie  activity.  The  operational  DOD  ionosphere  specification 
model  uses  both  boundaries  as  key  inputs.  Another  well  known  example  of  auroral  oval 
boundary  utility  is  the  high  latitude  radar  surveillance  systems  where  the  auroral  boundary 
location  is  used  to  help  discriminate  radar  clutter  from  images  of  projectile  objects.  HF  and 
radio  systems  near  the  auroral  zone  are  greatly  affected  by  auroral  precipitation  which  causes 
enhanced  ionization  in  the  ionospheric  E-region.  Both  the  equatorward  and  poleward 
boundaries  of  the  auroral  zone  are  operationally  derived  from  particle  flux  measurements 
acquired  by  DMSP  satellites. 

The  Hardy  Auroral  electron  flux  model  has  been  very  useful  for  modeling  high  latitude 
auroral  oval  boundary  since  its  initial  development  in  1985  by  Hardy  et  al.  (1985)  using  DMSP 
particle  flux  data.  The  original  Hardy  model  referred  to  as  H-85  here  is  extensively  used  for 
studying  auroral  precipitation  phenomena  in  the  space  weather  community.  It  was  upgraded 
in  2008  by  expanding  the  amount  of  DMSP  particle  spectrometer  data  used  in  the  analysis 
(Hardy  et  ah,  2008).  The  upgraded  Hardy  model,  referred  here  as  H-08,  turns  out  to  be  a 
large  database  and  lacks  functional  representation.  As  a  consequence  no  utility  has  been 
developed  for  it  yet.  To  explore  its  utilities,  we  attempted  to  derive  a  functional 
representation.  In  our  first  attempt  we  reduced  median  auroral  fluxes  from  the  H-08  model, 
and  conducted  a  regression  analysis  to  fit  the  median  auroral  fluxes  using  real  spherical 
harmonics  functions.  The  reduced  model  of  median  auroral  fluxes  we  obtained  is  named  as  H- 
08MV  to  distinguish  it  from  the  earlier  models.  The  purpose  of  this  report  is  to  describe  our 
preliminary  results.  We  first  review  the  H-85  and  H-08  models  in  the  next  section.  Next  in 
Section  3  we  compare  auroral  oval  boundaries  as  determined  by  these  two  Hardy  models. 
Section  4  describes  the  regression  analysis  using  real  spherical  harmonics.  In  the  Appendix  we 
list  tables  of  expansion  coefficients  in  the  functional  representation.  Finally  concluding 
remarks  are  given  in  Section  5. 

The  Hardy  Auroral  electron  flux  model  (H-85)  currently  used  in  the  space  weather 
community  was  originally  derived  from  precipitating  particle  measurements  made  by  the 
Special  Sensor  for  Precipitating  Particles,  version  4  (SSJ4)  electrostatic  analyzers  flown  on  the 
F6  and  F7  satellites  of  the  Defense  Meteorological  Satellite  Program  (DMSP)  (Hardy  et  al, 
1985;  Hardy  et  al,  1987).  A  total  of  2  years  of  ion  and  electron  data  (1983  data  for  F6  and 
1984  data  for  F7)  were  used  to  develop  the  H-85.  The  SSJ/4  analyzers  measured  electron  and 
ion  spectra  in  the  local  satellite  zenith  direction  once  per  second  over  20  channels  spanning  an 
energy  range  from  30  eV  to  30  keV.  Statistical  hemispheric  particle  precipitation  maps  were 
created  for  a  range  of  different  magnetospheric  activity  levels  defined  by  the  Kp  index.  The 
particle  precipitation  map  is  defined  in  the  coordinate  of  corrected  geomagnetic  latitude  (MFat) 
and  magnetic  local  time  (MET).  The  high-latitude  region  grid  was  defined  by  30  zones  in 
MFat  between  50“  and  90“  and  48  half-hour  zones  in  MET.  Each  latitude  zone  was  2“  wide 
between  50°  and  60“  and  1°  wide  between  60“  and  80°  latitude.  Measurements  in  each  of  the 
energy  channels  were  accumulated  to  build  a  one-second  average  of  the  differential  number 
flux  spectra  for  each  bin.  The  resulting  average  spectra  were  then  extrapolated  to  100  keV  and 
integrated  over  energy  to  yield  integral  number  flux  and  integral  energy  flux.  The  differential 
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number  flux  spectra  were  assumed  to  be  isotropic  in  the  integration.  The  average  energy  was 
calculated  as  ratio  of  integral  energy  flux  to  integral  number  flux.  Figure  1  produced  by  the 
AF-Geospace  Software  (Hilmer,  2010)  illustrates  the  precipitating  electron  energy  flux  map  in 
the  Northern  hemisphere  for  Kp  =  3. 


Figure  1.  An  example  of  the  precipitating  electron  energy  flux  map  in  the 
Northern  hemisphere  for  Kp=3  produced  by  AF-Geospace  software. 

The  green  line  specifies  the  auroral  oval  boundary. 

The  H-85  model  was  upgraded  to  FI-08  in  2008  by  re-analyzing  historical  DMSP  particle 
spectrometer  data  from  1983  to  2005  from  a  series  of  9  DMSP  satellites  (Hardy  et  al.,  2008). 
The  statistical  study  was  conducted  using  more  than  600  million  energetic  electron  spectra 
measured  by  the  SSJ4  sensor  to  obtain  total  number  fluxes  Jtot  in  the  unit  of  #/cm^-s-sr  and 
energy  fluxes  JEtot  in  the  unit  ofkeV/cm^-s-sr.  As  a  big  database  the  H-08  model  is  difficult  to 
use.  It  does  not  have  explicit  algorithm  for  calculating  the  auroral  oval  boundary.  It  divides 
the  ranges  of  total  number  flux  and  average  energy  in  logarithmic  scale  into  26X26  probability 
maps  for  each  bin  of  Kp,  MLat,  and  MET.  The  H-08  model  has  Kp  integer  values  from  0  to  9, 
MEat  from  43  to  89  and  MET  from  0  to  23. 

Some  improvements  over  the  H-85  model  are; 

(1)  H-08  divides  data  into  bins  1°  by  1  hr  in  magnetic  latitude  and  magnetic  local  time, 
respectively.  The  H-08  latitude  bin  is  1°  from  geomagnetic  latitude  50°  to  the  pole,  whereas 
the  H-85  latitude  bin  is  1°  wide  only  between  60°  and  80°.  However  the  original  H-85  model 
has  a  magnetic  local  time  bin  width  V2  hour  smaller  than  H-08. 

(2)  H-08  provides  probabilities  of  encountering  different  levels  of  total  number  flux  hot,  total 
energy  flux  JEtot  and  average  energy  (Eave  =  JEtot/  Jtot),  whereas  no  probability  statistics  are 
available  from  H-85.  H-08  demonstrates  that  within  all  MLat-MET  bins  probabilities  for 
encountering  different  levels  of  Jtot,  JEtot,  and  Eave  are  lognormal  distributed.  Often  more  than 
one  population  can  be  found. 
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(3)  Both  the  H-85  and  H-08  models  represent  data  according  to  the  level  of  the  Kp  index. 
The  H-85  electron  flux  characteristics  were  determined  for  each  whole  number  value  of  Kp 
from  0  to  5  and  for  Kp  >  6-.  The  H-08  model  specifies  the  electron  flux  probability  for  the 
whole  Kp  range  from  0  to  9.  However  the  statistics  are  poor  for  Kp  >  8  because  of  its  small 
sampling  size. 

(4)  The  methodologies  used  in  the  two  Hardy  models  are  different.  The  H-85  model  was 
developed  by  first  taking  the  average  number  flux  spectra  in  each  energy  bin  and  then 
integrating  them  to  obtain  Jtot  and  JEtot.  In  contrast  Hardy  et  al.  (2008)  integrated  first  number 
and  energy  fluxes  over  energy  and  then  performed  statistics  to  derive  their  probability. 

Figure  2  illustrates  a  representative  output  of  the  H-08  model  for  the  evening  (20  MET) 
sector  at  auroral  (72°,  69°,  and  65°)  and  subauroral  (60°)  magnetic  latitudes.  Contours  of 
constant  probability  densities  are  plotted  as  functions  of  log  JEtot  (abscissa)  and  log  Eave 
(ordinate),  where  log  is  the  10-base  logarithmic.  With  these  coordinates  a  single  lognormal 
distribution  would  appear  as  a  circular  region  with  an  intensity  peak  at  the  center.  As  indicated 
in  Figure  2,  probability  distributions  at  auroral  latitudes  were  represented  by  multiple 
lognormal  distributions  in  different  ranges  of  Eave  and  JEtot.  The  top  panel  of  Figure  2  (MEat  = 
72°)  shows  clear  evidence  of  two  populations  that  are  distinct  in  both  Eave  and  JEtot.  Peaks 
occur  in  Eave  near  0.2  and  1  keV.  In  the  subauroral  example  (bottom  panel,  60°),  Eave- JEtot 
characteristics  are  those  of  conjugate  photoelectrons  at  Eave  <  100  eV  that  were  encountered 
during  -'20%  of  the  DMSP  orbits.  The  third  panel  (65°)  also  shows  clear  conjugate 
photoelectron  population  at  the  lower  left  corner. 

From  the  occurrence  probability  maps  such  as  those  shown  in  the  first  column  of  Figure  2, 
the  H-08  model  can  be  used  to  yield  the  average-energy-fiux  probability  PF(log  JEtot)  that 
specifies  the  probability  of  detecting  average  auroral  energy  flux  at  the  level  of  JEtot  in  an 
MEat-MET  and  Kp  bin.  The  average-energy-fiux  probability  PF(log  JEtot)  is  obtained  by 
integrating  the  model  lognormal  distributions  over  Eave-  The  second  column  in  Figure  2 
illustrates  PF(log  JEtot)  as  a  function  of  log  JEtot  corresponding  to  the  first  column.  The 
average-energy-flux  probability  PF(log  JEtot)  is  comprised  of  several  components  as  plotted  in 
colored  lines.  Purple,  yellow  and  green  traces  indicate  fitted  distributions  of  the  primary, 
secondary  and  tertiary  populations,  respectively.  Blue  lines  are  sums  of  fitted  probability 
distributions.  Red  curves  represent  measured  probability  distributions.  The  top  and  third 
panels  indicate  that  PE  (log  JEtot)  is  represented  by  bimodal  lognormal  probability  distributions. 
The  second  and  fourth  panels  on  the  other  hand  suggest  PF(log  JEtot)  is  dominated  by  a  single 
component.  Similarly  the  H-08  model  can  be  used  to  specify  the  probability  for  detecting 
average  energy  at  a  specified  level  of  Eave  (Figure  2,  right  column). 
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Probability  Density  Functions 


Figure  2.  Four  representative  samples  of  the  probability  distributions  output  from 
the  Hardy  H-08  model  along  the  MLT  =  20  meridian  for  Kp  =  3. 

Their  magnetic  latitudes  are  shown  at  the  left  side.  Left  column:  contours  of 
constant  probability  densities  as  functions  of  log  JEtot  (abscissa)  and  log  Eave 
(ordinate).  Middle  column:  average  energy  flux  probability  PF(log  JEtot)  as  a 
function  of  log  JEtot.  Right  column:  average-energy  probability  PF(log  Eave)  as  a 
function  of  log  Eave-  (from  Efardy  et  al.  2008).  Red  curves  represent  measured 
probability  distributions.  Purple,  yellow  and  green  traces  indicate  fitted 
distributions  of  the  primary,  secondary  and  tertiary  populations,  respectively. 
Blue  lines  are  sums  of fitted  probability  distributions. 

2.  Auroral  Oval  Boundary  Model 

The  H-85  model  can  provide  a  map  of  auroral  boundary  at  all  longitudes  for  a  given 
geomagnetic  activity.  It  contains  an  algorithm  to  specify  the  equatorward  auroral  oval 
boundary.  The  auroral  oval  boundary  algorithm  is  not  based  on  the  spatial  map  of  model 
electron  fluxes.  Instead  it  was  deduced  from  a  separate  statistics  study  [Gussenhoven  et  al., 
1981;  Gussenhoven  et  al,  1983].  The  equatorward  boundary  was  determined  for  each  DMSP 
pass  from  latitudinal  profiles  of  the  electron  integral  number  flux.  A  large  number  of  DMSP 
boundary  crossings  (over  200,000  passes)  in  the  interval  1983-1990  were  included  in  the 
statistics  [Madden  and  Gussenhoven  1990].  Because  the  auroral  energy  flux  and  boundary 
models  are  based  on  two  independent  statistical  studies,  the  equatorward  edge  of  the  H-85 
energy  flux  distribution  is  close  to  but  do  not  always  coincide  with  the  auroral  oval  boundary 
determined  in  the  model.  For  example  the  green  line  in  Figure  1  that  specifies  the  auroral  oval 
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boundary  from  the  model  lies  outside  the  area  of  eleetron  preeipitating  fluxes.  These 
boundaries  would  be  eloser  with  a  different  ehoiee  of  threshold  flux  value. 

The  Gussenhoven  et  al.  [6]  regression  analysis  linearly  fit  the  equatorward  boundary 
versus  the  Kp  index  for  24  hourly  seetors  in  MLT.  The  equatorward  boundary  A  was 
expressed  in  a  linear  equation  of  Kp, 

A  =  Aq  +  aKp  (1) 

where  Ao  and  a  are  the  intereept  and  slope,  respeetively,  for  eaeh  magnetie  loeal  time  bin. 
Figure  3  illustrates  the  variation  of  the  auroral  equatorward  boundary  A  versus  Kp  in  the  2300- 
2400  hour  bin.  This  figure  shows  that  A  deereases  from  about  67°  when  Kp  =  1  to  less  than 
55°  when  Kp  =  7. 


Figure  3.  Mean  values  (solid  line)  and  standard  deviation  (vertieal  bar)  in 
eaeh  Kp  bin  as  a  function  of  Kp  for  the  2300-2400  MLT  sector  (from 
Gussenhoven  et  a/.  1983). 


Recently  the  equatorward  boundary  statistics  have  been  updated  in  the  AF-Geospace 
software  with  the  extensive  DMSP  auroral  particle  data.  Table  Al  in  Appendix  A  lists  the 
updated  values  of  Ao  and  a  for  24  hourly  sectors.  Note  that  Table  Al  contains  regression 
coefficients  for  a  few  MLT  bins  that  were  originally  unavailable  in  the  original  work  of 
Gussenhoven  et  al.  (1983). 

Based  on  the  H-85  auroral  boundary  model.  Madden  and  Gussenhoven  (1990)  developed 
an  auroral  boundary  index  as  the  equivalent  midnight  boundary  after  removing  local  time 
variation.  The  temporal  variation  of  the  midnight  auroral  boundary  index  follows  auroral 
activity  and  can  be  used  to  predict  global  changes  in  the  auroral  oval.  It  was  demonstrated  that 
the  midnight  auroral  boundary  index  more  accurately  reflects  magnetic  activity  than  the  Kp 
index  at  time  scale  much  shorter  than  the  Kp’s  three-hour  increment  (Madden  and 
Gussenhoven,  1990). 
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As  a  big  database  the  H-08  model  is  diffieult  to  use.  It  does  not  have  an  explieit  algorithm 
for  ealeulating  the  auroral  oval  boundary.  It  divides  the  ranges  of  total  number  flux  and 
average  energy  in  logarithmie  seale  into  26X26  probability  maps  for  eaeh  bin  of  Kp,  MLat,  and 
MLT.  The  H-08  model  has  Kp  integer  values  from  0  to  9,  MLat  from  43  to  89  and  MLT  from  0 
to  23. 

Similar  to  Gussenhoven  et  al.  (1983),  we  manually  seanned  the  H-08  probability  maps  to 
determine  the  threshold  value  of  MLat  where  the  oeeurrenee  probability  of  deteeting  eleetron 
flux  drops  rapidly  to  zero.  Speoifieally  the  two  eriteria  for  determining  an  auroral  boundary 
were  electron  number  flux  greater  than  10  #/cm  -sec-sr  and  the  encounter  probability  greater 
than  0.5.  The  low  and  high  latitudinal  boundaries  were  determined  for  each  MLT  and  each  Kp 
from  0  to  6.  Due  to  insufficient  coverage  for  Kp  >  6  the  probability  maps  do  not  have  good 
statistics  to  yield  a  good  average  model  of  auroral  boundary  at  Kp  >  6.  The  magnetic  latitudes 
of  the  auroral  boundaries  are  fitted  by  sinusoidal  functions  of  MLT  for  each  Kp  bin  as 


Z7I— 4  .  . 

(  nnx  ,  .  nnx\ 

and 

Zn=3  V  , 

(  nnx  ,  j  .  nnx\ 


(2) 

(3) 


where  and  Ah  are  the  MLat  of  the  low  and  high  latitude  boundary,  respectively,  and  x  is 
MLT  expressed  in  degrees  (x  =  15*MLT).  The  Al  expansion  coefficients  and  b„  up  to  n  = 
4  harmonics  are  given  in  Appendix  B  (Table  Bl)  for  Kp  from  0  to  6.  Because  the  high 
latitude  boundary  Ah  appears  to  vary  gradually  with  MLT,  Ah  is  fitted  only  up  to  n  =  3 
harmonics.  Its  coefficients  c„  and  dn  are  listed  in  Appendix  B  (Table  B2).  From  these  tables 
the  auroral  boundary  MLat  for  fractional  Kp  value  can  be  specified  by  linearly  interpolation. 

Figure  4  compares  the  H-08  low  latitude  auroral  boundary  Al  determined  from  Equation 
(2)  with  the  H-85  equatorward  boundary  A  expressed  in  Equation  (1)  for  Kp  =  1,3,  and  5. 
This  figure  shows  that  Al  (solid  line)  is  slightly  lower  than  A  (red  triangle)  in  the  1000-1600 
MET  hours.  Eor  other  METs,  Al  and  A  are  about  the  same. 
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Figure  4.  Comparison  of  auroral  oval  boundaries  between  the  Fl-85  and  H-08 
models. 

The  solid  line  in  each  Kp  panel  represents  the  low  latitude  auroral  boundary  Al 
determined  from  the  H-08  probability  maps.  The  H-85  equatorward  auroral 
boundary  from  Equation  (1)  is  plotted  as  solid  triangle. 

3.  H-08  Median  Values  (H-08MV)  Model 

An  auroral  specification  capability  based  on  the  H-08  model  would  have  a  wealth  of  statistical 
information  that  could  potentially  provide  new  and  improved  information.  For  instance  more 
accurate  confidence  levels  could  be  provided  to  system  operators  and  better  error  estimates 
could  be  passed  to  models  assimilating  auroral  fluxes.  Here  we  discuss  a  first  attempt  to  make 
the  H-08  more  accessible;  this  model  was  studied  in  the  validation  effort  described  in 
Section  5. 

Since  the  H-08  model  is  a  large  database  and  lacks  functional  representation,  no  utility  has 
been  developed  for  it  yet.  To  explore  its  utilities,  we  attempted  to  derive  a  functional 
representation  of  the  H-08  model.  We  first  determined  the  medians  of  PF(log  JEtoO  and  PF(log 
Gave)  for  each  MLat-MLT-Kp  bin  from  their  H-08  model  occurrence  probability  maps, 
respectively.  Figure  5  illustrates  the  obtained  median  energy  flux  log  JEtot  as  a  function  of 
MEat  at  MET  =  0  for  Kp  from  0  to  6. 
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Figure  5.  Median  energy  flux  from  the  FI-08MV  model  as  a 
function  of  magnetic  latitude  at  MLT  =  0  for  Kp  from  0  to  6. 

In  order  to  expand  log  JEtot  and  log  Eave  in  terms  of  spherical  harmonics  functions,  we 
define  a  normalized  variable  0  between  the  lower  and  upper  auroral  boundaries  {Al  and  Aif)  as 

_  n{MLat-Ao) 

where  Ai  is  defined  as  {Al  +  Ah)I2.  Note  that  0  varies  from  -nil  to  nil  as  MEat  varies  from  Al 
to  Ah.  Similarly  we  define  (()  =  7i:MLT/12  that  varies  from  0  to  2n  to  replace  the  MLT  angle 
from  midnight.  Eor  each  Kp  multivariate  regression  analysis  was  performed  to  fit  median 
values  of  JEtot  and  Eave  versus  0  and  ((),  respectively,  by  real  spherical  harmonics  yik{Q,  (j)) 
(www.wikipedia.org)  as: 

(f)  =  Ci^yikid.  <f)  (5) 

and 

Eavei.O.(p)  =  (6) 

(j)  is  MLT  expressed  as  angle  from  midnight  in  radian,  Cik  and  Dik  are  the  expansion  coefficients 
determined  from  regression  analysis  using  up  to  4th  order  real  spherical  harmonics  ym  functions 
(/  from  0  to  4,  and  k  from  0  to  21).  Tables  Cl  and  C2  in  Appendix  C  list  the  regression 
coefficients  Cik  and  Dik,  respectively,  for  Kp  =  0  to  6.  The  function  representation  of 
Equations  (5)  and  (6)  along  with  lookup  tables  of  Cik  and  Dit  constitute  the  reduced  model  of 
median  auroral  fluxes  and  energy,  which  we  named  it  the  H-08  Median  Value  (H-08MV) 
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model  to  distinguish  it  from  the  earlier  Hardy  models.  The  H-08MV  representation  allows 
straightforward  ealeulation  of  median  auroral  energy  flux  and  median  energy  at  given  MLat 
and  integer  Kp.  Interpolation  is  applied  for  fractional  Kp.  A  Fortran  computer  code  for  the 
H-08MV  model  is  available  from  the  authors  upon  request. 

Figure  6  shows  contour  plots  of  the  median  auroral  energy  flux  and  median  energy  for  Kp 
=  0  to  6.  The  H-08MV  contours  indicate  that  the  low  latitude  boundary  of  auroral  electron 
precipitation  extends  from  70°  to  55°  as  Kp  increases  from  0  to  6.  Both  JEtot  and  Eave  medians 
increase  with  Kp. 


K  p  :  0  K  p  :  1  K  p  :  2  K  p  =  3  K  p  ;  4  K  p  =  5  K  p  =  6 


□ 
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Eigure  6.  Contour  plots  of  the  median  auroral  energy  flux  and  median 
energy  from  the  H-08MV  model  for  Kp  =  0  to  6. 

4.  Discussion 

The  H-85  and  H-08MV  models  have  been  recently  examined  in  the  validation  study  of  auroral 
oval  models  by  Eane  (2012)  in  which  several  available  models  specifying  the  auroral  oval 
boundary  were  validated  with  DMSP  satellite  measurements.  Eane  (2012)  represented  the 
equatorward  boundary  as  the  location  at  which  energy  flux  measured  by  DMSP  satellites 
exceeded  a  fixed  threshold  of  0.4  erg/cm2/s.  The  observed  auroral  boundary  locations  were 
compared  with  predictions  from  H-85,  H-08MV,  Ovation  Prime  model,  Space  Weather 
Modeling  Framework  model,  and  Assimilative  Mapping  of  Ionospheric  Electrodynamics 
(AMIE)  model. 

The  validation  study  used  data  obtained  from  3  DMSP  satellite  passes  during  a  total  of  28 
days  in  the  period  of  2000  to  2008.  These  dates  included  13  days  covering  five  prominent 
geomagnetic  storms  and  15  days  with  low  and  moderate  Kp  values.  A  total  of  6,297  DMSP 
satellite  passes  were  included  in  this  study.  Unfortunately  due  to  limited  data  availability  no 
validation  study  was  performed  for  critical  local  times  from  22  to  03  MET. 

Lane  (2012)  studied  the  model  prediction  efficiency  (PE)  by  evaluating  the  coefficient  of 
efficiency  defined  as 


PE  =  1- 


Zf=i(yi-y)" 


(7) 
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where  N  is  the  number  of  data  points,  yt  is  the  i-th  data  point,  y  is  the  mean  of  the  data  set,  and 
Xi  is  the  model  prediction  value  correspond  to  j,.  The  PE  score  measures  the  percentage  of 
variance  in  the  observed  DMSP  data  captured  by  the  model.  Thus,  a  PE  score  of  1  indicates 
perfect  prediction  by  the  model.  An  efficiency  of  0  (PE  =  0)  indicates  that  the  model 
prediction  is  as  accurate  as  the  mean  of  the  observed  data. 

Eane  (2012)  found  the  Ovation  Prime  model  to  have  the  highest  overall  prediction 
efficiency  score  (PE  =  0.55).  The  PE  score  is  0.51  for  H-85  and  0.45  for  H-08MV,  respectively. 
These  PE  scores  imply  that  the  empirical  auroral  models  have  considerable  variance  in  their 
predictions.  This  conclusion  is  supported  by  the  earlier  work  of  Gussenhoven  et  al.  (1983), 
which  showed  that  the  Hardy  auroral  boundary  model  has  a  standard  deviation  of  about  one 
degree  for  Kp  <  3  and  larger  values  for  higher  Kp  (see  Eigure  3). 

5.  Conclusion 

Auroral  oval  boundary  has  been  considered  a  critical  requirement  for  DMSP  mission  and  space 
situational  awareness.  Auroral  particle  precipitation  fluxes  are  essential  environment  data 
records  acquired  by  DMSP  satellites  to  meet  the  requirement.  Both  the  equatorward  and 
poleward  boundaries  of  the  auroral  zone  are  the  derived  EDRs  from  particle  flux  data.  The 
threshold  requirement  for  auroral  boundary  is  the  identification  of  the  auroral  boundaries  along 
the  satellite  path.  Although  both  the  equatorward  and  poleward  boundaries  are  of  interest,  the 
equatorward  boundary  is  of  greater  priority  than  the  poleward  boundary.  Specification  of  the 
poleward  boundary  is  usually  considered  as  an  objective  requirement  since  its  determination  is 
often  ambiguous. 

The  original  Hardy  auroral  oval  H-85  model  has  been  extensively  used  for  modeling  high 
latitude  phenomena  and  space  weather  operation  since  its  initial  development  in  1985.  The  H- 
85  model  gives  the  average  spectra  of  particle  fluxes  from  which  mean  values  of  various 
quantities  like  energy  flux  and  mean  energy  are  derived.  One  important  objective  of  the  model 
is  to  map  the  auroral  boundary  at  all  longitude  for  a  given  geomagnetic  activity.  The  H-85 
model  has  achieved  a  high  level  of  maturity  with  explicit  algorithms  to  predict  auroral 
boundary  location. 

The  H-85  model  has  recently  been  extended  into  a  probability  model  H-08  with  a  much 
larger  collection  of  new  DMSP  observations  since  1985.  Its  improvement  mitigates  the  H-85 
deficiency  by  providing  variance  information.  As  a  large  database  the  upgraded  Hardy  H-08 
model  is  nevertheless  difficult  to  use  and  requires  a  functional  representation.  The  present 
study  has  demonstrated  a  feasible  approach  to  represent  the  H-08  model  using  real  spherical 
harmonics  functions.  Briefly  the  functional  representation  we  explored  so  far  models  only 
median  energy  flux  and  median  energy.  To  expand  utility  of  the  H-08  model  we  are  currently 
extending  the  present  study  to  represent  other  model  features.  Eurther  improvement  will 
greatly  enhance  its  usefulness  for  space  weather  prediction. 
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Appendix  A.  H-85  Model  Regression  Coefficients  for  the  Equatorward  Auroral  Boundary 

Table  A-1 .  H-85  model  regression  eoeffieients  for  the  equatorward  auroral  boundary  A 


K  =  Ao  +  aKp 


MLT 

Number* 

Ac 

a 

cc' 

0000-0100 

17235 

66.9 

-1.79 

-0.76 

0100-0200 

7234 

66.7 

-1.57 

-0.70 

0200-0300 

8282 

66.6 

-1.64 

-0.75 

0300-0400 

16153 

67.0 

-1.81 

-0.78 

0400-0500 

35610 

67.0 

-1.88 

-0.78 

0500-0600 

67812 

67.7 

-1.81 

-0.78 

0600-0700 

72618 

68.2 

-1.71 

-0.73 

0700-0800 

56758 

68.8 

-1.58 

-0.69 

0800-0900 

65537 

69.7 

-1.46 

-0.65 

0900-1000 

61515 

70.0 

-1.19 

-0.53 

1000-1100 

31844 

70.2 

-0.85 

-0.41 

1100-1200 

14340 

70.2 

-0.41 

-0.22 

1200-1300 

7954 

71.3 

-0.21 

-0.12 

1300-1400 

7390 

73.1 

-0.44 

-0.23 

1400-1500 

8313 

72.8 

-0.62 

-0.34 

1500-1600 

11733 

73.4 

-1.09 

-0.59 

1600-1700 

26411 

73.5 

-1.38 

-0.73 

1700-1800 

45591 

72.6 

-1.51 

-0.76 

1800-1900 

61810 

71.7 

-1.66 

-0.80 

1900-2000 

66432 

71.0 

-1.83 

-0.82 

2000-2100 

79626 

70.0 

-1.81 

-0.81 

2100-2200 

63606 

69.2 

-1.94 

-0.82 

2200-2300 

44349 

68.1 

-1.92 

-0.83 

2300-2400 

32577 

67.2 

-1.83 

-0.81 

*Number  of  boundaries  used  in  eaeh  regression 
tCorrelation  eoeffieient  between  A  and  Kp. 
Souree:  AF-Geospaee  software  (Hilmer,  2010). 
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Appendix  B.  Coefficients  for  the  H-08  Model  Boundary 


Table  B-1.  Coefficients  for  the  H-08  model  low  latitude  boundary  Al  in  Equation  (2) 


Index 

Kp=0 

Kp=l 

Kp=2 

Kp=3 

Kp=4 

Kp=5 

Kp=6 

ao 

1.1758 

1.1504 

1.1271 

1.1068 

1.0800 

1.0544 

1.0378 

ai 

-0.0222 

-0.0484 

-0.0554 

-0.0744 

-0.0851 

-0.0800 

-0.0915 

bi 

-0.0114 

-0.0255 

-0.0383 

-0.0388 

-0.0275 

-0.0216 

-0.0176 

^2 

-0.0114 

-0.0204 

-0.0145 

-0.0097 

-0.0020 

0.0046 

0.0130 

hj 

-0.0067 

-0.0080 

0.0115 

0.0011 

0.0119 

0.0175 

0.0129 

a3 

-0.0011 

0.0023 

0.0083 

0.0049 

0.0095 

0.0094 

-0.0013 

b3 

0.0010 

0.0009 

-0.0013 

-0.0020 

-0.0077 

-0.0039 

-0.0049 

a4 

0.0006 

-0.0002 

-0.0059 

0.0007 

-0.0037 

-0.0067 

-0.0013 

b4 

-0.0013 

-0.0027 

0.0027 

0.0013 

0.0010 

0.0038 

0.0050 

Table  B-2.  Coefficients  for  the  H-08  high  latitude  boundary  Ah  in  Equation  (3) 


Index 

Kp=0 

Kp=l 

Kp=2 

Kp=3 

Kp=4 

Kp=5 

Kp=6 

Co 

1.3860 

1.3969 

1.3890 

1.3898 

1.3911 

1.3796 

1.3867 

Cl 

-0.0829 

-0.0538 

-0.0398 

-0.0174 

-0.0018 

0.0150 

0.0125 

di 

0.0116 

0.0222 

0.0070 

0.0140 

0.0195 

0.0315 

0.0183 

C2 

-0.0169 

-0.0026 

-0.0100 

-0.0011 

-0.0020 

-0.0041 

-0.0095 

^2 

0.0051 

-0.0038 

-0.0027 

0.0077 

0.0104 

0.0085 

0.0217 

C3 

0.0039 

0.0060 

-0.0009 

-0.0009 

-0.0015 

-0.0095 

-0.0076 

d3 

-0.0022 

-0.0107 

0.0024 

0.0012 

-0.0098 

0.0022 

0.0024 
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Appendix  C.  Regression  Coefficients  for  the  H-08MV  Model 


Table  C-1.  Regression  eoefficients  Cik  in  Equation  (5) 


Index 

Kp=0 

Kp=l 

Kp=2 

Kp=3 

Kp=4 

Kp=5 

Kp=6 

O 

o 

o 

2.5933e+01 

2.6419e+01 

2.7018e+01 

2.7331e+01 

2.7642e+01 

2.7752e+01 

2.7679e+01 

Cio 

4.2699e-02 

2.9665e-02 

9.6545e-02 

2.3478e-01 

3.9974e-01 

5.9882e-01 

8.1152e-01 

C,1 

2.9979e-01 

2.7641e-01 

4.1384e-01 

4.2845e-01 

3.7436e-01 

2.9723e-01 

3.7766e-01 

C,2 

-4.027  le-01 

-6.1018e-01 

-7.7352e-01 

-7.9939e-01 

-6.7944e-01 

-5.0470e-01 

-7.0573e-01 

C20 

-4.7319e-01 

-6.3424e-01 

-8.3750e-01 

-9.8990e-01 

-1.0881e+00 

-1.2260e+00 

-1.2178e+00 

C21 

-8.2667e-02 

-4.6545e-02 

1.7062e-02 

-2.9322e-02 

-1.5787e-01 

-2.6405e-01 

-1.1977e-01 

C22 

-9.6017e-02 

-1.0683e-01 

-1.1731e-01 

-2.3017e-02 

1.0577e-01 

-2.5746e-02 

-1.3074e-01 

C23 

-1.0602e-01 

-1.7478e-01 

-2.6656e-01 

-4.6299e-01 

-5.3215e-01 

-5.0706e-01 

-5.3579e-01 

C24 

4.6688e-02 

-7.0096e-02 

2.5682e-03 

-7.5863e-02 

-7.7337e-02 

-2.0728e-01 

-2.8664e-01 

C30 

4.1655e-01 

4.8183e-01 

5.0093e-01 

5.0559e-01 

5.0510e-01 

4.6267e-01, 

4.4003e-01, 

C31 

7.8009e-02 

1.0694e-01 

6.5848e-02 

7.7903e-02 

1.2989e-01 

1.6437e-01 

8.4770e-02 

C32 

8.7785e-02 

9.901  le-02 

1.3824e-01 

1.5905e-01 

1.2522e-01 

5.091  le-02 

3.5080e-02 

C33 

-1.5315e-01 

-1.4413e-01 

-7.2309e-02 

-7.5693e-03 

8.6032e-02 

-1.2168e-02 

1.7777e-01 

C34 

-8.7330e-03 

2.9115e-01 

2.0476e-01 

2.6477e-01 

1.5019e-01 

2.5081e-01 

2.4226e-01 

C35 

-2.5404e-01 

-2.8757e-01 

-3.0073e-01 

-1.9864e-01 

-1.3104e-01 

-1.5396e-01 

-1.0487e-01 

C36 

2.6598e-02 

-1.1223e-01 

-2.0551e-01 

-4.0589e-01 

-3.2158e-01 

-2.0497e-01 

-1.6785e-01 

C40 

-7.0876e-02 

-6.5799e-02 

4.5761e-02 

-5.5574e-02 

-4.9123e-02 

-9.5350e-02 

-1.4024e-01 

C41 

3.7265e-02 

1.2984e-01 

1.4247e-01 

2.2153e-01 

3.0379e-01 

1.6304e-01 

1.4159e-01 

C42 

1.7147e-01 

3.0712e-01 

2.5850e-01 

3.1306e-01 

2.1320e-01 

2.5894e-01 

2.1720e-01 

C43 

3.4809e-02 

-6.3848e-02 

5.2888e-03 

5.0161e-02 

2.1075e-01 

1.6739e-01 

1.5175e-01 

C44 

1.6392e-02 

-1.4506e-01 

-3.8450e-03 

7.3798e-03 

-6.4012e-02 

-3.6353e-03 

-1.8067e-02 

C45 

4.4232e-02 

6.4640e-03 

9.8167e-02 

2.0145e-01 

-2.6796e-04 

1.2838e-01 

-4.6818e-03 

C46 

5.6995e-03 

4.498  le-02 

-7.7958e-02 

-1.9563e-01 

-1.7600e-01 

-2.1337e-01 

-2.3249e-01 

C47 

4.9388e-02 

2.8978e-02 

8.2980e-02 

5.1840e-02 

2.7555e-02 

5.6134e-02 

5.5444e-02 

C48 

-3.7796e-02 

-1.0546e-02 

-6.2586e-03 

-4.8818e-02 

-4.8359e-02 

-7.6932e-03 

-3.4973e-02 

Note:  Coeffieient  Cik  matehes  the  sequenee  with  real  spherieal  harmonics  yik  in  Appendix  D. 
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Table  C-2.  Regression  eoefficients  Dik  in  Equation  (6) 


Index 

Kp=0 

Kp=l 

Kp=2 

Kp=3 

Kp=4 

Kp=5 

Kp=6 

O 

o 

p 

-1.7300e+00 

-1.4493e+00 

-1.1307E+00 

-8.9516E-01 

-8.2792E-01 

-8.8079E-01 

-1.0317E+00 

Dio 

4.8954e-02 

1.9877e-01 

4.8066E-01 

7.4207E-01 

8.3337E-01 

9.8394E-01 

1.0719E+00 

Di, 

3.4805e-01 

4.5990e-01 

5.9400E-01 

5.8077E-01 

4.9866E-01 

3.8492E-01 

3.1449E-01 

Di2 

-4.7046e-01 

-5.9127e-01 

-5.1359E-01 

-5.0540E-01 

-4.6081E-01 

-3.0305E-01 

-3.4897E-01 

D20 

-1.0936e+00 

-1.1012e+00 

-1.1802E+00 

-1.1446E+00 

-1.2044E+00 

-1.1934E+00 

-9.7272E-01 

D21 

-1.7796e-01 

-2.8930e-01 

-1.1049E-01 

-1.9860E-01 

-2.8294E-01 

-4.0946E-01 

-2.5643E-01 

D22 

2.4059e-02 

7.9987e-02 

-4.4199E-02 

-7.7953E-02 

-3.9396E-02 

-4.6639E-02 

-1.8697E-01 

D23 

-9.4771e-02 

-2.5964e-01 

-3.7743E-01 

-5.2293E-01 

-5.0814E-01 

-4.6820E-01 

-3.7312E-01 

D24 

-2.3983e-02 

-6.5280e-02 

-1.4453E-01 

-2.5212E-01 

-3.0463E-01 

-3.5296E-01 

-3.5772E-01 

D30 

5.2075e-01 

6.4961e-01 

6.2316E-01 

5.9178E-01 

5.7904E-01 

4.8779E-01 

4.3183E-01 

D3, 

1.6378e-02 

2.6462e-01 

2.7206E-01 

2.7523E-01 

3.0326E-01 

2.7434E-01 

4.0121E-02 

D32 

-7.0786e-02 

-8.4989e-02 

1.5623E-01 

2.0750E-01 

1.9372E-01 

2.0480E-01 

1.2278E-01 

D33 

7.1226e-02 

,-3.4012e-02 

1.0322E-01 

1.3190E-01 

1.0930E-01 

2.5162E-02 

1.4691E-01 

D34 

8.2478e-02 

4.2146e-01 

2.5173E-01 

3.4261E-01 

4.6852E-01 

4.2942E-01 

4.5759E-01 

D35 

2.2487e-02 

-6.4025e-03 

-1.9770E-01 

-1.4019E-01 

-8.5956E-02 

-9.2246E-02 

-9.775  lE-03 

D36 

2.6611e-01 

9.3448e-02 

1.7421E-01 

1.0094E-01 

-6.5296E-03 

-2.0879E-02 

-9.0104E-02 

D40 

2.4341e-01 

1.4175e-01 

1.4508E-01 

4.3554E-02 

5.3648E-02 

3.2140E-02 

-4.8405E-02 

D41 

-5.0369e-02 

4.7840e-02 

4.5021E-02 

2.9833E-02 

7.5888E-04 

-8.3912E-02 

-4.6096E-02 

D42 

4.7256e-02 

1.5081e-01 

1.1735E-01 

1.6831E-01 

1.4838E-01 

2.1080E-01 

1.2729E-01 

D43 

-7.6145e-02 

-3.0821e-02 

-4.4963E-02 

-2.2661E-02 

1.6409E-03 

8.4459E-02 

3.7660E-03 

D44 

4.6876e-02 

-7.3577e-03 

1.8354E-01 

1.8351E-01 

1.4177E-01 

1.6503E-01 

2.2308E-02 

D45 

1.2508e-02 

6.2900e-02 

-1.1056E-02 

2.1242E-02 

-1.6333E-01 

-5.7537E-02 

-2.2938E-02 

D46 

3.9841e-02 

-2.1440e-02 

-1.1262E-01 

-1.6892E-01 

-8.9952E-02 

-1.5480E-01 

-1.5214E-01 

D47 

6.2767e-03 

8.3762e-02 

1.3513E-02 

-5.9936E-02 

-1.0670E-01 

-9.8622E-02 

-8.2133E-02 

D48 

2.0683e-02 

-6.3000e-02 

-1.6166E-01 

-1.4319E-01 

-1.4308E-01 

-8.8456E-02 

5.9848E-03 

Note;  Coefficient  Dik  matches  the  sequence  with  real  spherical  harmonics  ym  in  Appendix  D. 


Approved  for  public  release;  distribution  is  unlimited. 

16 


Appendix  D.  Formulas  of  Real  Spherical  Harmonics 


The  formulas  of  the  real  spherieal  harmonics  7^(0,  (j))  used  in  this  study  are  based  on  “Table  of 
spherical  harmonics”  in  www.wikipedia.org.  These  formulas  that  are  related  to  the  complex 
spherical  harmonics  7/^(0,  (}))  are  expressed  in  the  Cartesian  coordinate  (x,  y,  z,  and  r).  We 
transform  0  and  (()  to  the  Cartesian  coordinate  through  the  usual  spherical-to-Cartesian 
coordinate  transformation  with  r  =  1 : 


X  —  rsinO  coscj) 
y  =  r  sin  9  sincj) 
z  =  r  cos  6 


Real  spherical  harmonics  with  /  =  0  is: 


Too 


Real  spherical  harmonics  with  /  =  1  are: 


Tio 


Til 


Real  spherical  harmonics  with  1  =  2  are: 


720 


721  =  77  — -r? 


722  -  n 


723  -  n  I  ■  „9 


724 


X 

r 

}Att 

r 

h  ■ 

7 

i47r 

r 

z 

|47r 

r 

1^- 

— X' 

Jtt 

15 

yz 

r 

y,2 

15 

zx 

h 

y2 

h 

xy 

‘Z 

y,2 

■ +  2z^ 
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Real  spherical  harmonics  with  /  =  3  are: 


1  7  z(2z^  —  3x^  —  3y^) 


y3i 


733  -  4 


734-2 


■j^3 

[35 

(3x2 

-7^)7 

1 27r 

■j^3 

I35 

(x2  - 

-  3y2)x 

jzrr 

■j^3 

105 

(x2 

-y^)z 

71 

■j^3 

|l05 

xyz 

71 

■j^3 

1  21  y(4z^— x^— y^) 


736 


1  21  x(4z^— x^— y^) 
4,  27r 
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Real  spherical  harmonics  with  /  =  4  are: 


3  1  (35z^  -  30z^H  -  3r^) 


740  —  ^(1  _  ■  ■ 

16^  n 


741 


3  5  xz{lz^  —  3r‘^) 


2  _  Q^2>, 


4j27r 


y‘i 

2  _ 


3  5  yz(7z^  —  3r^) 


y'-i 


3  5  (x"^  —  y"^)(7z"^  —  3r"^) 


y^ 

2  „2> 


3  5  xyilz^  —r^') 


744  —  .  ■ 

4.  TT 


2  _  Q-,,2>, 


745 


746 


3  35  (x"^  —  3y"^)xz 
4^271 

3  (3x^  —  y^)yz 

AAZtc 


3  35  x^(x^  —  3y^)  —  y^(3x^  —  y^) 


747  —  T7  “  ■ 
16  ^  TC 


3  35  xy(x"^— y"^) 


748  —  7 

4^  TT 
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